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<D • Abstract 



We introduce a new method to propagate uncertainties in the beam shapes used to measure the cosmic microwave background to 
cosmological parameters determined from those measurements. The method, called Markov Chain Beam Randomization (MCBR), 
randomly samples from a set of templates or functions that describe the beam uncertainties. The method is much faster than direct 
numerical integration over systematic 'nuisance' parameters, and is not restricted to simple, idealized cases as is analytic marginaliza- 
tion. It does not assume the data are normally distributed, and does not require Gaussian priors on the specific systematic uncertainties. 
f \ ' We show that MCBR properly accounts for and provides the marginalized errors of the parameters. The method can be generalized and 

used to propagate any systematic uncertainties for which a set of templates is available. We apply the method to the Planck satellite, 
and consider future experiments. Beam measurement errors should have a small effect on cosmological parameters as long as the 
CIh' beam fitting is performed after removal of 1 // noise. 
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, ^3- 1. Introduction 



Observations of the cosmic microwave background (CMB) can be interpreted only in light of a detailed knowledge of the angular 
response of the instrument to radiation, i.e., the shapes of the 'beams.' It is almost always the case that the beams from single- 
■ aperture telescopes (but not interferometers) can be approximated as two-dimensional Gaussians. It is never the case that a gaussian 
I/-) \ approximation provides an adequate description of the beams of an experiment that measures the CMB with high signal-to-noise 
i ratio. If the beams were known perfectly, their effects on the data could be calculated perfectly, if painfully. Unfortunately, beams 
' are never known perfectly, and among the outstanding issues for any CMB experiment are how to optimize the beams in the first 
| place, and how to control and account for beam uncertainties in the data analysis. 

The effects of beam uncertainties can be analyzed in maps, power spectra, and cosmological parameters determined from the 
' data. Each has benefits. Because cosmological parameters are a key product of any experiment, and because they are sensitive to 
f^i , extremely small effects impossible to detect pixel by pixel, they are particularly valuable. Historically, however, calculation of the 
r* ■ effects of beam uncertainties on cosmological parameters has been done either analytically, which requires over-simplified beam 
shapes, or numerically, at great computational cost. 

We introduce in this paper a method for calculating the effects of beam uncertainties on cosmological parameters determined 
^ ■ from CMB observations that is both fast and flexible. It requires only that beam uncertainties, or for that matter any other systematic 
| effect, can be represented by a set of functions or templates, which could be obtained from Monte Carlo simulations. It does not 
assume that the data themselves are Gaussian-distributed, or that the uncertainties have Gaussian priors. 

In § II we describe the method, called Markov Chain Beam Randomization or MCBR, and we show that the MCBR tech- 
nique produces correct marginalized errors. §111 summarizes the beam fitting procedure developed in a previous paper 



(N 
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( Huffenberger et al. 20091. §IV describes the implementation of MCBR. In § VI we apply the method to the Planck experiment, 



and consider future experiments. 



2. MCBR: Markov Chain Beam Randomization 



In the past, marginalization over systematic parameters has been carried out either numerically or analytically ( IBridle et al. 2002t ; 
both methods are currently implemented in cosmomc dLewis and B ridle 2002). Assuming likelihoods are Gaussian one typically 
has a marginalization of the form: 



Loc daP(a)exp[-(a\ - d) r N _1 (av - d)/2] 



(1) 
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where d is the theoretical (predicted) data minus the observed data and a\ is an approximate template describing the correction 
applied to the predicted data due to systematics, N is the noise covariance matrix, and P(a) is the prior. The marginalization is given 
by: 

-21nL = d r 1ST 1 =— = d + ln^N^'v) + c (2) 



where c is a constant. If v is independent of the data and parameters then L oc e ~ x ?ff^ 2 , with: 



x - eff = d T \N- 1 - yrN _ ly j d = xtest-fu (3) 

In the case of beam uncertainties, the analytic approach is feasible only if the beams are assumed to be Gaussian. This is not 
realistic. 

It is customary to characterize anisotropies in the Cosmic Microwave Background by their angular power spectrum, Cf for both 
temperature and polarization. Cf is a 3 x 3 matrix for T (temperature) and E or B (grad-type or curl-type polarization): 



C,= 



pTT pTE pTB \ 

L L L 

pTE pEE pEB 

p V. p p 

pTB Aeb Abb 

*-t W W 



(4) 



Hereafter, for the sake of simplicity, most equations will refer to the angular power spectrum, Q, for a single component, say 
temperature. The telescope beam smooths the anisotropies, supressing power at higher multipoles. We refer to the ratio of the 
measured power spectrum of the sky and our true power spectrum as the transfer function, £( = B 2 . Here we assume the beam 
transfer functions are the same for temperature and polarization. 

To obtain unbiased estimates of the parameters that characterize the cosmology, we must repair this suppression based on 
knowledge of the beam. Uncertainties in the beam propagate into uncertainties in the cosmological parameters. 

We assume that the beam uncertainties can be described by a set of functions or templates, taken here to be the set of transfer 
functions obtained by the beam fitting procedure described in § [3] These templates are given in multipole space by: 

SJ = {B' ( f = (B e x r t ) 2 (5) 

where the ratios rc represent the possible deviations from the true fiducial beam. We choose the beam transfer function randomly 
from the set of N simulations (here N = 1280) for each step of the Markov Chain Monte Carlo when probing the cosmological 
parameters space. This means that at each step of the chain the theoretical power spectrum, Cf, is multiplied by the randomly chosen 
beam, S r ( . We assume all transfer functions in the set are equally probable. 

To estimate constraints on cosmological parameters, we need to compare the model with the data via a chosen Likelihood 
and an algorithm to sample cosmological parameters. Here we make use of the package cosmomc. To incorporate MCBR we mod- 
ify cosmomc to enable the usage of a random S r ( for each theoretical model generated with CAMB dLewis et al. 20001) or PICO 
dFendt and Wandeir2 006). This is done by modifying the cmbdata module of the cosmomc code. 

We start by creating simulated datasets with noise properties specific to the instrument under consideration, in our case Planck 
and an example of a future experiment (see § @). These simulated datasets are given in terms of the angular Power Spectrum C° bs : 

r obs = C wmaP Sf + Nf (6) 

where Nc is the noise power spectrum and Be = B 2 is the beam transfer function and C™ ap is the A CDM spectrum best-fitting 

current WMAP data. In the case of a symmetric Gaussian beam, B(x) — — ^ exp J^j-J, so that B( = e However, the Planck 

beams are not adequately represented by Gaussians. Instead, we use realistic beams calculated from a full diffraction analysis of the 
telescope using GRASP9 dSandri et al. 20021 ISandri et al 20091 iMaffei et al. 20091 1 Yurchenko et al. 20041 . 

As our purpose here is to introduce and validate the MCBR method it suffices to assume full-sky coverage. Considerations of 
realistic complications (such as cut-sky, foregrounds, etc.) is deferred to a future publication. Our purpose here is to establish 
the relative importance of propagating beam errors to cosmological parameters rather than to make comprehensive predictions 

C { = Cf s ) dBond, Jafle and Knox 2000) : 



for Planck. Hereafter to compare the observed dataset, Cl bs , with theoretical models we use the exact full-sky likelihood (with 



- 2 In L(C e \C t ) = (K + 1) (in |Q| + Tr (c, C^ 1 )) , (7) 

i.e., the Inverse Wishart distribution for Temperature and Polarization. In cosmomc this distribution is coded in function ChiSqExact 
dLewis 20051 ). We analyse these datasets with a modified version of this function, built to include the MCBR procedure in the code. 
The C( of the theoretical model is given by: 

C e = C ( x S r e (8) 

where T> r ( . is the randomly chosen transfer function. To incorporate both the beam and the uniform white noise in the likelihood 
expression one should replace: 
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C e -> Cf s (9) 
Q -> C' e h xS' ( + N[ (10) 

where C° ( bs is given by Equation[6] Cf is the theoretical power spectrum computed e.g. by CAMB, and B r ( is the randomly chosen 
beam transfer function. 

In the MCBR scheme, sampling of the beam templates is equivalent to sampling from the proposal distribution. The Metropolis- 
Hastings algorithm accepts the move from 9„ to 9 n+ \ in the Markov chain by evaluating the ratio: 

P(Qn+l)q(Qn+\,Qn) ,^y. 
P(0„)q(0„,0 n+1 ) 

where P is the posterior distribution we wish to sample from and q is the proposal distribution. We draw the proposal at position 
6„ of the parameter space from q(6„+\,6„). Here 6 = (6 cp ,6b), with 6 cp the subset of cosmological parameters and 9b the beam 
parameter. The joint proposal density for 9 factors into 

q(fin+l,Qn) - qcp{6cp,n+l,6cp,n)qb(6b,n+l,Qb,n), (12) 

where C p, n +\ refers to (9 cp ) n+ \ and &b, n +\ to {9b)„+\. Now, we take the qb{0b, n +\,Qb,n) to be the posterior distribution of the beam 
parameters given the beam fitting data (in our case the Jupiter beam fitting data (see § [3}), i.e., 

qb{0b,n+u0b,n) = Pi(%„|beamdata) (13) 

Furthermore 

P(9 n+l ) = P e ^„ + i|mapdata)/m,H+i|beamdata) (14) 

(for instance in our study here P cp (#|mapdata) = L(mapdata|0 f; ,, 9b)p cp (9 cp ) where L is the Likelihood given in Equation [7] and p cp 
the prior on cosmological parameters.) Hence the ratio in Equation [TTlbecomes : 

Pcp(@n+l)qcp(@cp,n+l> @cp,n) (15) 
Pcp(0n)Qcp\ficp,ni @cp,n+\) 

as Pb and qb cancel out. 

Hence random sampling from the set of beam templates at each step of the Markov chain is equivalent to sampling from a 
proposal density that, by construction, is identical to the posterior distribution of the beam parameters given the beam fitting data. 

To illustrate how the MCBR procedure works, we give here the steps followed in our analysis (see section [5]). We start by com- 
paring two cases: 

1. cosmomc run with the 'true' fiducial beam transfer alone, St. 

2. cosmomc run with the MCBR procedure for the set of beam transfer functions, S' r obtained from the beam fitting step. 
To this end: 

- We generate a simulated data set using as fiducial the 'true' beam transfer Sf 

- We analyse this simulated data set with cosmomc, including in the code just the effect of the S{ (ie the theoretical Q is multiplied 
by the true beam transfer, Sf ) 

- We analyse this simulated data set with a modified version of cosmomc in which the theoretical C7 is multiplied by the randomly 
chosen beam, S' ( at each step of the chain — i.e., with in-built MCBR 

The MCBR method can be used to propagate any systematic uncertainties that can be characterized by a set of templates. We 
turn these into multiplicative and additive corrections to the C(, encode the corrected Cf into the likelihood, and randomly sample 
from the set of templates at the Markov Chain Monte Carlo step of parameter estimation. The data do not have to be normally 
distributed. Furthermore, unlike analytic marginalization, the method does not require Gaussian priors on the uncertainties. 



2.1. Validation 

To demonstrate that that the MCBR technique gives correct marginalized errors, we compared the results given by MCBR to those from 
a 'brute force' cosmomc calculation in which the beam was taken as another parameter. We did this for three simulated datasets, the 
first generated using the 'true' beam transfer function S( , the second and third using beam transfer functions that were chosen to be 
mildly and extremely far from the true one, respectively. We simplified the test cases by assuming that the beam was a symmetric 
7' (FWHM) Gaussian, with 4% variations of the fwhm of the beam. 

The brute force calculation was done by probing the beam parameter space in cosmomc in the same way as for any of the other 
parameters, and considering the default proposal density already implemented in cosmomc. The beam parameter is included by 
transforming the theoretical C((9 cp ) output by CAMB at each Markov chain step into 

Ce(0c P , Ob) = C c (9 cp ) x B" + N u (16) 
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Figure 1. Marginalized parameter constraints for Planck 143 GHz with 7 beam with 4% variations, for the analysis with the 'true' 
reference fiducial beam using beamparameter approach (black) and MCBR (red), the blue line is the analysis o f the same dataset 
without including the beam uncertainty. 



where fwhm,, is the width of B", the Gaussian beam currently sampled. This theoretical C((9 cp , 6b) is used in the Likelihood expres- 
sion. To move the Markov chain to the next position in parameter space we use the default proposal density in cosmomc, usually 
an Af-d Gaussian. The proposed new point is accepted/rejected following the same prescription used for the other cosmological 
parameters. A final marginalized distribution of the beam is output along with the other cosmological parameter constraints. 

For the MCBR calculation we analysed the simulated data with a modified version of cosmomc in which the theoretical Q is 
multiplied by the randomly chosen beam transfer, 3 r ( at each step of the chain. 

The results are plotted in Figures Q] |2j and[3] In all cases, we find same parameter distributions for both methods. As expected, 
the extreme deviated beam results in a biased estimation of parameters, especially n s , but equally for both the 'beamparameter' and 
the MCBR procedures. 



3. Beam fits and transfer function ensembles 



We characterize the beam uncertainty for Planck with a Monte Carlo ensemble of transfer functions ( Huffenberger et al. 2009 ) 
generated by repeated simulation of Jupiter observations using the detector noise and pointing errors expected before flight. Each 
realization yields a representative transfer function. The beams are calculated with GRASP9 dSandri et al. 20021 ISandri et al 2 009 
IMaffei et al. 20 09 Yurchenk o et al. 20041 1. and we employ two methods of beam reconstruction to reproduce them from the planet 
scans. The first uses a rigid linearized parametric model; the second expands the beam in orthogonal functions (see lRocha et al. 20011 
for a previous application of such functions in CMB analysis). Figure [5] shows the nominal Gaussian beams with blue-book fwhm 
values to that of the fiducial realistic Grasp beams based on a Gaussian fit (see tabled). From the beam reconstruction procedure 



presented in Huffenberger et al. 2009 we obtain the ratio of the power spectrum as corrected with the fitted beam to the power 
spectrum as it should have been corrected by the true beam. In Figure|4]we display lines which bound 68% of the ensemble transfer 
functions for Planck channels. 

The simulation of repeated Jupiter calibrations id done in such way that each template is an unbiased estimator of the true 
template. But in real life, they could be a biased estimator (for instance the Planck pointing error could bias the beam function 
always in the same direction). This prompted us to consider the runs presented in Section l5T2l 

The Beam fitting is applied to data with white +1// noise, and to destriped data, i.e., after application e.g of a "de- 
striping" mapmaking code which removes almost all of the effects of 1 // noise IPoutanen et al. 20061 lAshdown et al. 2007al 
lAshdown et al. 2007bl lAshdown et al. 20091 . We use realistic Grasp beams and the parametric model of the reconstructed beams 
(the results with non-parametric model will be presented in a future paper). Figure|6]shows extreme and mild beam transfer functions 
for the Planck 70 GHz, 100 GHz, 143 GHz and 217 GHz channels obtained from the beam fitting procedure applied to destriped data 
(hence containing a a very low level of 1 // residuals). For comparison purposes we plot in Figure [7] these functions obtained from 
data with a white and 1 // noise background. We also plot in Figures[8]the normalized histograms of the ratios, r?, for singe multi- 
poles 6 = 500, 1000, 1500, 2000. For all channels the distributions are slightly skewed and get broader with increasing multipole I 
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Figure 2. Marginalized parameter constraints for Planck 143 GHz with 7 beam with 4% variations, for the analysis with the 'true' 
reference fiducial beam (black) and for the mildly deviated beam transfer, SJ - (B{ x r"" ld ) 2 , using beamparameter approach 
(blue) and MCBR (red), both beamparameter and MCBR give same distributions 




Figure 3. Marginalized parameter constraints for Planck 143 GHz with 7 beam with 4% variations, for the analysis with the 'true' 
reference fiducial beam (black) and for the extremely deviated fiducial beam, & ( xt - (B{ x r e ( xt ) 2 , using beamparameter approach 
(blue) and MCBR (red), both beamparameter and MCBR give same distributions. 
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Figure 4. At each multipole, 68% of the fitted Monte Carlo transfer functions recover spectra closer to the true power spectrum 
than the indicated line. Left: parametric model. Right: non-parametric model based on orthogonal functions, where the flexibility 
requires less knowledge of the beam, but yields larger errors. 
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Figure 5. Nominal Gaussian blue-book beams (dotted line) vs Fiducial realistic Grasp9 beams based on a Gaussian fit (solid line) 
for70GHz (black), 100GHz (red), 143GHz (green) and 217GHz(blue). 



for each channel. We can also compare the probability of the mildly deviated and extremely deviated transfer functions used in § [5] 
For instance for 70 GHz for I - 1000 the mild function is 10% probable while the extreme function is approximately 100 times 
less likely. The maximum variation for transfer function ratios, r 2 c , is of the order 2% for 70 GHz (0.5% for 100GHz) for destriped 
data, while for white and l/f noise data with no attempt at destriping it increases to ~ 30% for the 70 GHz channel (~ 2.5% for 
100 GHz). 

In § |5]we infer that the parameter constraints from beams obtained with destriped data are slightly worse but very close to those 
obtained with a white noise background as expected. 



4. Analysis: from Beam transfer function uncertainties to parameter estimation 

To propagate the beam measurement errors to parameters we apply the MCBR method following the procedure described in § [2] We 
make use of the beam transfer functions obtained with the beam fitting described in § [3] For this purpose we use a modified version 
of cosmomc with built-in MCBR step as described in § [3J We consider a set of five chains. The convergence diagnostic is based on 
the Gelman and Rubin statistic, as usual in the field. Following MCBR, we choose randomly the beam transfer function (from the 
set of 1280 simulations) for each step of the Markov Chain. We sample a six-dimensional set of cosmological parameters, with flat 
priors: the physical baryon and Cold Dark Matter densities, u>b = Q/,/r and to c = Q c /j 2 ; the ratio of the sound horizon to the angular 
diameter distance at decoupling, 6 S ; the scalar spectral index its; the overall normalization of the spectrum /og[10 10 A] at k = 0.05 
Mpc 1 (hereafter As), and the optical depth to reionization t. We use a cosmic age top-hat prior 10Gyr< to < 20Gyr, consider 
purely adiabatic initial conditions only, we impose flatness, and we treat the dark energy component as a cosmological constant. 

We create simulated datasets with the noise properties of the Planck 70, 100, 143 and 217 GHz (IPlanck B lue Book 2005) chan- 
nels, as well as one example of a future experiment. For the latter we considered the noise levels of Epic 150 GHz ( IBock et al. 20081 1. 
We take as our cosmological model the best fit of WMAP lyr: Q b h z = 0.02238; Q c /i 2 = 0.11061; H = 71.992; r = 0.110267; 
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Figure 6. Extreme (left) and Mild (right) beam transfer functions for the Planck 70 GHz (black), 100 GHz (red), 143 GHz (green) 
and 217 GHz (blue) channels obtained from beam fitting applied to destriped data. 
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Figure 7. Extreme (left) and Mild (right) beam transfer functions for the Planck 70 GHz (black), 100 GHz (red), 143 GHz (green) 
and 217 GHz (blue) channels obtained from beam fitting applied to data with white +1// noise. 




Figure 8. Normalized distributions of the beam transfer functions, S ( = (B' ( ) 2 for multipoles I = 500, 1000, 1500, 2000 for 70 GHz 
(black), 100GHz (red), 143 GHz (green), and 217GHz(blue) obtained from the beam fitting on destriped data. 



«s = 0.95820; and As = 3.0824 (Spergel et al. 2003). These simulated datasets are given in terms of the angular Power Spectrum 
C°^ s as described in § [2] We compute the noise Ne = (AT x fwhm) 2 for Planck and Epic from the sensitivity AT /T and the nominal 
fwhm of the beam assuming a Gaussian profile (tabulated in Table[TJ. In Figure|9]we plot the theoretical model vs. the noise levels 
for each channel considered. Results from this analysis are given in § [5] 
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Experiment Channel FWHM AT/T FWHM (grasp) 

Planck 70 W 47 IT 

100 10' 2.5 9.22' 

143 7.1' 2.2 6.49' 

217 5.0' 4.8 4.48' 

Epic-CS 150 5.0' 0.81 

Table 1. Planck (Planck Blue Boo k 2005 ) and Epic (Boc k et al. 2008t experimental specifications. Channel frequency is given in 
GHz, FWHM in arcminutes and noise in 10~ 6 . The last column gives the fwhm of our fiducial beams based on a Gaussian fit to the 
realistic GRASP beams. 
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Figure 9. CMB angular power spectrum (best fit of WMAP lyr, black line) and noise levels for Planck: 70 Gflz(black), 
100 Gtfzfred), 143 G//z(green), 217 Gtfz(blue) and for Epic 150 G#z(cian). 



5. Results 

5. 1 . Results: effect of beam uncertainties 

Figures \W\ QT| [12] Q/3] show the marginalized parameter constraints for Planck in three cases: without beam uncertainty (i.e., 
considering the true fiducial beam) (black); with beam uncertainty using the beam transfer functions obtained using the destriped 
data (red); and in the presence of 1 // noise (blue). Table [2] gives the input cosmological parameters and the mean values and 
marginalized 68% confidence limits obtained after accounting for the beam errors. To facilitate comparisons, Figure[14]shows these 
same marginalized constraints for n s and A s , the parameters where the largest differences are seen between the three cases. 

Equivalent results for a more sensitive polarization experiment — Epic 150 GHz — are plotted in Figure Q3J corresponding pa- 
rameter values are given in table [2] 

The most noticeable effect of beam uncertainties is to widen the marginal distributions of some parameters, especially n s , for 
uncertainties obtained in the presence of 1 // noise but without destriping. A s , and Q/,/? 2 are also affected. In this case the distribution 
of the fitted beam transfer functions is wider than that obtained from white noise or destriped data as shown in Figures [6] and [7] 

Define <t c /,, the width of the distribution when beam errors are marginalised by applying MCBR, and <rw, the width of the distri- 
bution for the simulated data convolved with the fiducial beam (with no beam errors included). Figure [16] shows the enhancement 
factor, cr c h/cr re f, for parameters n s and A s for beams fitted on data with white + l/f noise. For example, at 100 GHz the distributions 
of n s and A s widen by 25% and 11%, respectively. 

This widening is much reduced by the use of destriping techniques. For example, with destriping the uncertainties in the beams 
are 0.5% for 100 GHz at 6 = 1500, which translates into an increase of parameter uncertainties of 0.1%. Without destriping, the 
uncertainties on A s at 70GHz and on n s at 100GHz increase by 21% and 25%, respectively, for beams fitted on white + l/f noise 
data. This is a convincing demonstration of the relevance and power of destriping techniques in reducing the effect of 1 // noise for 
Planck. 

5.2. Results: effect of assuming a wrong fiducial beam 

To illustrate the effect of incorrect beam assumptions we calculated parameters assuming a mildly and then an extremely 'wrong' 
beam. Specifically, we generated three simulated datasets, using: Be ; a mildly wrong beam ( B'" dd ; and an extremely wrong beam S e ( xl . 
We analyzed these datasets with the modified version of cosmomc with MCBR built in, and compared the cosmological parameters 
from the run with for the 'true' fiducial beam !B( to those of both the mild and extreme deviated beams 

Figures [171181191201 show marginalized parameter constraints from 70GHz, 100GHz, 143 GHz, and 217 GHz, respectively, 
on destriped data. We see that assuming an extreme beam deviation in the simulated data results in a biased estimation of some 
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Channel 


Parameter 


no beam uncertainty 


destriped 


white noise+ 1 / f 


Planck 70 GHz 


Q b h 2 


0.22393 ± 0.00035 


n 994.01 


-+- n nnfns 


0.22394 ± 0.00036 




£l c h 2 


0.1106 ± 0.0027 


n 1 1 n^ 


+ n nfi97 


0.1106 ± 0.0029 




e 


1.0428 ±0.0010 




+ n nn 1 n 

It U.UulU 
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Table 2. Mean values and marginalized 68% c.l. limits using the fiducial beam: analysis without beam uncertainty (column 3), 
accounting the beam uncertainty from destriped data (column 4) and from the data with white and 1 // noise (column 5). 




Figure 10. Marginalized parameter constraints for Planck 70 GHz without beam uncertainty (black), marginalized over the beam 
uncertainty via MCBR considering the destriped data (red), and in the presence of white noise +1// noise (blue). 
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Figure 11. Marginalized parameter constraints for Planck 100 GHz without beam uncertainty (black), marginalized over the beam 
uncertainty via MCBR considering the destriped data(red) and in the presence of white noise +1 // noise (blue). 
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Figure 12. Marginalized parameter constraints for Planck 143 GHz without beam uncertainty (black), marginalized over the beam 
uncertainty MCBR considering the destriped data (red) and in the presence of white noise + l/f noise (blue). 
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Figure 13. Marginalized parameter constraints for Planck 217 GHz without beam uncertainty (black), marginalized over the beam 
uncertainty MCBR considering the destriped data (red) and in the presence of white noise + l/f noise (blue). 
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Figure 14. Marginalized constraints for the most impacted parameters, n s and A s , for Planck channels 70GHz, 100GHz and 143GHz, 
without beam uncertainty (black), marginalized over the beam uncertainty considering the destriped data (red) and in the presence 
of white noise +1// noise (blue). 
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Figure 15. Marginalized parameter constraints for a future experiment with Epic 150 GHz specifications without beam uncertainty 
(black), marginalized over the beam uncertainty considering the destriped data (red) and in the presence of white noise +1// noise 
(blue). 
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Figure 16. Enhancement factor, <r c h/cr re f for n s (solid line) and A s (dashed line), where cr c h is the width of the distribution when 
beam errors are marginalized over by applying MCBR; cr re f is the width of the distribution for the simulated data convolved with the 
fiducial beam (no beams errors included), for beams fitted on data with white and 1 // noise background. 



parameters, particularly n s . This is mostly due to incomplete marginalization, as we do not encompass an adequate distribution of 
deviations from the chosen fitted transfer function. 

For comparison, Figures |2T1 and l22l show marginalized parameter constraints for the 100GHz and 143GHz channels, respec- 
tively, on data that have not been destriped. 

Figure [23] shows the bias in n s and A s as a function of the extreme beams fitted on destriped data. We consider the error on B( 
given by fa - 1) for i = l/cr representing the sigma of the beam. The corresponding values are given in Table [3] For example for 
100GHz an uncertainty of the extreme beam transfer function b 2 ( for I — 810 of 0.1% bias the likelihood by 0.3<x and 0.13cr for 
n s and A s respectively. A beam transfer function known up to 0.02% will bias n s by 0. lo\ If we had not taken into account the beam 
uncertainties, then the same deviation in the transfer functions would have biased n s by as much as 0.4cr, as can be inferred from 
Figure[l6] The inadequacy of a likelihood that does not integrate the beam uncertainties is mentioned in ( Huffenberger et al. 2009 1. 
There a simplified analysis of noisier data (only 1 horn) with all parameters except n s fixed indicated that limiting the bias to O.lcr 
would require knowledge of b 2 c to 0.04% where it has fallen to 1% of peak (£ 1900 for 100 GHz). In our analysis here we see that 
at { 1900 an uncertainty of 0.5% for the extreme function would bias n s by 0.3cr, while a mild deviation of the order 0.2% would 
produce a bias below 0.05<x (see Table[2]). Hence a beam deviation five times that reported in ( |Huffenb erger et al. 2009 1 would bias 
n„ by less than O.lcr. This improvement is mostly due to properly marginalizing over the beam uncertainties via the MCBR method. 



ch I rj bias/<r(n s ) bias/cr(A s ) 

~1Q 579 1.00217 0.1890 0.2846 

100 810 1.000928 0.3282 0.1240 

143 1141 1.000982 0.3775 0.1055 

217 1620 1.00043 0.0929 0.0188 

Table 3. Bias on n s and A s in units of the error due to the deviation of the extreme function (r"') 2 at £ = l/cr, after MCBR, fitted on 
destriped data. For each Planck channel. 



6. Conclusions 

We have developed a fast new method, MCBR, to propagate beam uncertainties to parameter estimation. The method properly ac- 
counts for the marginalised errors in the parameters. A desirable feature of the method is that it makes minimal assumptions on 
beam uncertainties. For example, it does not assume the data are normally distributed, and, unlike other approaches such as analytic 
marginalization, it does not require Gaussian priors on the specific systematic uncertainty. Furthermore it accounts accurately for 
the shape of the beam as it makes use of beam uncertainty templates for such beams, hence there is no need for simplified a priori 
assumptions on their shapes. Finally MCBR can be generalized and used to propagate other systematic uncertainties, as long as a set 
of templates of such systematics is provided. 

From the study presented here on propagating the beam measurement errors to parameter estimation via the new MCBR method 
for Planck and for a future experiment, we conclude: 

- Removal of l/f noise residuals, by destriping or other techniques, is quite important. 

- The main impact of beam uncertainties is to widen the marginal distributions of some parameters (most notably n s ). 
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Figure 17. Marginalized parameter constraints for Planck 70 GHz with beam randomization MCBR: true beam (black), decreasing 
function for destriped data (red), increasing function for destriped data (blue ), mild deviation (solid line) and extreme deviation 
from the true beam (dotted line) 




Figure 18. Marginalized parameter constraints for Planck 100 GHz with beam randomization MCBR: true beam (black), decreasing 
function for destriped data (red), increasing function for destriped data (blue ), mild deviation (solid line) and extreme deviation 
from the true beam (dotted line) 
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Figure 19. Marginalized parameter constraints for Planck 143 GHz with beam randomization MCBR : true beam (black), decreasing 
function for destriped data (red), increasing function for destriped data (blue ), mild deviation (solid line) and extreme deviation 
from the true beam (dotted line) 
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Figure 20. Marginalized parameter constraints for Planck 217 GHz with beam randomization MCBR: true beam (black), decreasing 
function for destriped data (red), increasing function for destriped data (blue ), mild deviation (solid line) and extreme deviation 
from the true beam (dotted line) 
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Figure 21. Marginalized parameter constraints for Planck 100 GHz with beam randomization MCBR: true beam (black), decreasing 
function for white +l/f noise (red), increasing function for white +l/f noise (blue ), mild deviation (solid line) and extreme deviation 
from the true beam (dotted line) 
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Figure 22. Marginalized parameter constraints for Planck 143 GHz with beam randomization MCBR: true beam (black), decreasing 
function for white +l/f noise (red), increasing function for white +l/f noise (blue ), mild deviation (solid line) and extreme deviation 
from the true beam (dotted line) 
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Figure 23. Bias on « ,(solid line) and A. s (dotted line) in units of the error for the extreme beam functions, r"' for I = l/cr y after 
beam randomization MCBR, fitted on destriped data. For 70GHz (black), 100GHz (red), 143GHz (green) and 217GHz (blue). 



- Assuming as extreme beam deviation in the simulated data results in a biased estimation of some parameters (mainly of n s ) due 
to incomplete marginalization. 

- The parameters more noticeably impacted by beam uncertainties are: n s , Q/,/) 2 and A , 

These results demonstrate the relevance of applying destriping techniques on Planck data to remove 1/f noise. 

When the beam fitting is performed in destriped data the uncertainties on the beams for say 100GHz are at most of the order of 
0.5% for t = 1500 which translates into an increase of parameter uncertainties at most of the order of 0. 1 %. Instead the uncertainties 
on A s at 70GHz and on n s at 100GHz increases approximately by 20% and 26% respectively for beams fitted on white + 1 // noise 
data while it remains unaltered for white noise background alone. 

The effect of wrong assumptions on beam parameters will bias the parameter constraints only for extreme deviations from 
the true beam and hence for quite atypical circunstances. Considering the analysis performed on destriped data, at 100 GHz an 
uncertainty of the extreme beam transfer function at I = 810 of = 0.1% will bias the likelihood by 0.3cr and 0.13cr for n s and A s , 
respectively. A beam transfer function known to 0.02% will bias n s by 0. lcr. If we had not taken into account the beam uncertainties, 
then the same deviation in the transfer functions would have biased n s by as much as QAcr. To limit the bias in n s to less than O.lcr 
will require a knowledge of a mild deviated beam fr 2 to 0.2% where it has fallen to 1 percent. A mild deviated function gives rise to 
no observable bias (ie at most of the order 0.05<x). 

Therefore we expect only a small impact of beam measurement errors on cosmological parameter estimation as long as the beam 
fitting is performed on destriped data. 
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